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We obtain a query lower bound for quantum algorithms solving the phase estimation problem. 
Our analysis generalizes existing lower bound approaches to the case where the oracle Q is given by 
controlled powers Q p of Q, as it is for example in Shor's order finding algorithm. In this setting we 
will prove a f2(log V 6 ) lower bound for the number of applications of Q P1 , Q P2 , . . .. This bound is 
tight due to a matching upper bound. We obtain the lower bound using a new technique based on 
frequency analysis. 
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I. INTRODUCTION 

We study lower bounds for the phase estimation prob- 
lem. In this problem we are given a unitary transfor- 
mation Q as a black-box and we know that \q) is an 
eigenvector of Q, i.e. 



\q)=e M *\q), p6[0,l). 



(1) 



We want to determine the phase tp up to precision e. 

The quantum phase estimation algorithm approxi- 
mates tp given \q) and is the main building block in 
Shor's factoring algorithm, the counting algorithm, and 
the eigenvalue estimation algorithm [1-7]. 

The main element of this algorithm are controlled pow- 
ers of Q, which we define as follows. Let Q be a t qubit 
unitary transformation and \ip) an arbitrary t qubit state. 
For I = 1, . . . , c, and p £ N we define the c+i qubit trans- 
formation 



Wf{Q) \ Xl 



■x c ) \i>) = 



\xi..:x c ) \ip) 
\ Xl ...x c )QP\iP) 



xi =0 
x\ = 1 



(2) 



We call Wf(Q) a (controlled) power query. If the trans- 
formation Q is clear from the context, we will just write 



Wf = 



Wf(Q). 



This notation allows us to write the phase estimation 
algorithm in a compact form, see figure 1. The algorithm 
returns an approximation tp of the phase tp of \q). 
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FIG. 1: The quantum phase estimation algorithm in power 
query notation. H is the Hadamard gate, Wf a power query 
as in equation (2), and is the inverse quantum Fourier 
transform on T qubits. 



It is well known that T = C(loge _1 ) power queries 
suffice to approximate tp up to e. In this paper we study 
whether it is possible to improve on the performance of 
the phase estimation procedure, i.e., we ask what is the 
minimal number of applications of Wf to estimate tp up 
to e. 

Theorem 1. Any quantum algorithm estimating the 
phase tp of an eigenvector \q) of matrices Q up to pre- 
cision e, with Q from the class 

Q\ q ) t t = {Q '■ Q is a unitary t qubit transform, 

\q) is an eigenvector o/Q}. (3) 

has to use f2(logi) power queries. 

We prove Theorem 1 in section IV. 

The query cost of the phase estimation algorithm may 
be given by counting each application of Wf(Q) as p 
applications of Q, i.e. the query cost of this algorithm is 



1 



lT-1 



1. 



For certain problems, like order- finding, it is possible 
to exploit some knowledge about Q. Here Q\y) — 
\xy mod N) for a certain fixed x, and therefore 

Q v \y) = \x v y mod N) = \ (x 2 ^ 1 ^ 2 y mod N) 

is easy to compute by repeated squaring and modular 
multiplication, see e.g. [2]. In this case we can use power 

queries Wf with essentially the same cost as an applica- 
tion of Q. Thus we can execute the phase estimation al- 
gorithm with query cost T and have exponential speedup 
for the query cost: from 2 T — 1 to T. Let us stress that 
this speedup only applies if the cost for computing Q 2 ' 
is similar to that for computing Q. 



II. PRIOR WORK 
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Quantum query complexity has been important in 
quantum computing since Grover's search algorithm, 



2 



which is provably superior to classical algorithms [8, 9] 
in the number of queries. The first lower bound result 
was given in [10], which used an adversary argument. 

Our lower bound approach is based on the ideas of the 
"polynomial approach" of Beals et. al., [11-14]. Other 
approaches include the quantum adversary argument and 
its generalizations [15-17]. 

These approaches only cover problems concerning Boo- 
lean functions, so we have to extend them to numerical 
problems. This has been done through extensions of the 
polynomial method [18, 19] and has been applied widely 
to integration [20-22], path integration [23], approxima- 
tion [24-26], and eigenvalue estimation [6, 7]. 

In this paper we apply the approach of [19] to the phase 
estimation problem. Instead of using a maximum degree 
argument, which is not applicable in the case of arbitrary 
powers, we will develop a new lower bound technique 
based on frequency analysis. 



III. QUANTUM ALGORITHMS WITH POWER 
QUERIES FOR PHASE ESTIMATION 

We would like to derive lower bounds for any quantum 
algorithm with power queries that estimates the phase of 
a matrix Q for a given eigenvector \q). In other words 
the set of allowed inputs for our problem is 

Q\ q ) t t = {Q '■ Q i s a unitary t qubit transform, 

\q) is an eigenvector of Q}. 

We now give a framework that is general enough to 
allow us to analyze any algorithm with power queries 
WF(Q) that solves this problem. The most general algo- 
rithm will be of the following form: 

\l>W{Q)) = U T Wg(Q)U T -i...W£(Q)U \^ 0) ). (4) 

Here the Uo, U i, . . ., Ut are arbitrary but fixed c+t qubit 
unitary transformations and |?/;( )) a fixed c + t qubit 
state, for example |0) \q). In our analysis we neglect the 
cost to implement the Uj or to prepare |?/>( )). 

The varying parts of algorithm (4) arc the Wf 3 — 
Wf 3 {Q): power queries of Q for pj G N, lj = 1, . . . , c, 
and c G N arbitrary. A measurement of the final 
state \^ T HQ)) m the standard basis yields a state \k), 
k = 0, . . . , 2 c+t — 1, with probability Pk,Q, from which 
we get a solution <p(k) G [0, 1). If for all Q G Q\ q ).t the 
probability to get an e-estimate to the correct phase ip of 

Q 



E 

k:\\<p-<p(k)\\<e 



Pk,Q > 4, 



(5) 



then the algorithm (4) solves the phase estimation prob- 
lem to within e with probability | in T power queries. 

We are interested in the smallest number T such that 
a quantum power query algorithm of form (4) fulfills con- 
dition (5) for all Q G Q\ q ).t- 



IV. GENERAL CONTROLLED ARBITRARY 
POWER QUERIES 

We consider arbitrary powers pi, . . ., pt- This requires 
us to introduce a new proof technique. To illustrate the 
idea consider the phase estimation algorithm, performed 
as in figure 1 with c = T = 3 control qubits. 

(J--3 1 ® I)W*W 2 Wi(H® 3 ® J)|0) \q) . 

Let us trace through each of the steps in this algorithm 
(we neglect normalization factors). 

1. (|0) + |l) + |2) + |3> + ... + |7»|g) 

2. (|0> + c 2 ^ |1) + 1 2) + e 2 ^ |3) + . . . + e 2 ™* |7}) \q) 

3. (|0) + c 2 ^ |1) + e 2ni2 f |2> + . . . + e 2 ™ 3 ^ |7)) \q) 

The possible multiplicities j of ip in the coefficients 
e 27rw are 

J2 = {0,Pl,P2,Pl+P2} = {0,1,2,3}. 

4. (|0) + e 2 ^ |1> + e 2m2 ^ |2> + . . . + e 2m7 ^ |7)) \q) 
The possible multiplicities after this step are 

J3 = {j.i+P3:je^ 2 } = {0,l,...,7}. 

The final step, the inverse Fourier transform, does not 
depend on (p. It also does not change the possible mul- 
tiplicities of ip, but just creates linear combinations of 
them. Consider, e.g., the coefficient of the state |2): 

7 7 
^2 e -2Ki2j/a e 2xij<p | 2 ) |g) = ^ e 2irij(v.-l/4) ^) \ q ) ; 

which gives the probability pi (ip) of measuring |2): 



P2(<P) 



J=0 



7 

3,1=0 



,27r»C7'-0(v-V4) 



which is plotted in figure 2. Figure 2 shows that the 
probability that \2) is measured is high if ip is close to 
0.25, which is the value represented by |2). 

The figure indicates that the width of this probability 
peak depends on the frequencies present in the proba- 
bility function: higher frequencies allow sharper peaks. 
The goal of this paper is to prove that every halving of 
the width of the probability peak requires one additional 
step of the algorithm. 

The proof consists of three steps. The first is to 
quantify the influence of each additional application of 
Wf on the frequencies present in the probability func- 
tion (Theorem 2). Now consider a probability function 
as in figure 2. It must have a high peak > 3/4 with 
small width e and should be close to zero everywhere 
else. We will show that such a function requires a large 
range of frequencies to be present (lemma 3). Finally we 
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FIG. 2: (Color online) The probability of measuring the state 
1 2) depending on cp for the algorithm depicted in figure 1 with 
T = 2 and T = 3. 



will show an upper bound on the number of frequencies 
that a power query algorithm can generate with T power 
queries, which proves Theorem 1. 

In Theorem 2 consider a general Q, which acts on t 
qubits and therefore has 2* eigenvectors \\jj s ) with eigen- 
values e 27 ™ Vs . We assume that the eigenvectors \tp s ) are 
fixed, but that the eigenvalues change. We will prove 
that after T steps only coefficients like 

Q ,g27ri(j'iV3i+...+i 2 t»3 2 t) 

will occur. Here (Ji, . . . , J2*) is from the set Jt defined 
by the recursion 

J T+ l ■■= {(ill — ,32*), (il +PT+1, -J2*), 

(ii,-,i2* +pt+i) ■ (ii,-,i2*) e J T ) (6) 

and Jo = {(0,...,0)}. 

Theorem 2. Let Q be a unitary operation with eigen- 
vectors \ip s ) and corresponding eigenvalues e 27!ltp3 , s = 
1,...,2*. Let the \tp s ) be fixed and vary the phases 
tp s G [0, I). Any quantum algorithm with power queries 
Wf — Wf(Q), fixed unitary transformations Uj and 
starting state |^>(°)), can be written as 

U T W[ T T ...U 1 W^U \^)=J2 S P(^-^2 t )\k) (7) 



(T) 

for all ip s G [0, I), where the Sj, (<fii, ■ ■ ■ , </?2') are 
trigonometric polynomials of the following form: 



S { P(ip u ...,ip 2 t) 



E 



(T) 



a 2Tri(j 1 tp 1 + ...+j 2 tip 2 t) 



(8) 



with a k ^ - t ^ G C and Jt defined by (6). 



To shorten our proofs we will use the following short- 
hand notations. We say that a trigonometric polynomial 
S^\ip±, . . . , ifi2* ) is of powers Jr, indicated by the 
superscript, if it can be written as a sum over J7r as in (8). 
Furthermore we abbreviate the vectors j = {ji, ■ ■ ■ , J2* ) > 
= (fx, . • • , (^2* ) : an d define the following notation: 

j ■(f=jiip 1 + ... +j 2 *(f2 t - 

Proof of Theorem 2: To simplify the proof we choose a 
different basis instead of the standard basis. We divide 
the quantum state into a control part m) and an eigen- 
vector part \ip s ) and show that we can write 

2 C — 1 2* 

UtW^.-.W^Uo]^) = EE^($K^>> ( 9 ) 

m— s — 1 

for trigonometric polynomials of powers Jt-> 



(10) 



jeJ T 



The proof is by induction on the number of queries T. 
For T — we have no dependence on and Sm^sip) is 
just a constant since Q was not applied: 



U 



2 C -1 2 ^ 



^ (0) ) = E E^'^N^) i m '^> 



m— s—1 

2 C -1 2* 

m— s—1 

2 C -1 2* 

: EEE 

rn=0 g=l J £jb 



f (0) 

Sn,s,(0,...,0) 



m,s,j 



m,ip s ) 



Let now T be arbitrary and let equations (9) and (10) 
hold. If we apply Wf^^ to (9), only those states |m, ip s ) 
are affected for which the Ix+i-th. control bit is set, i.e. 
mi T+1 = 1- For these states we get 

|m, V.) = \m) Q Pr+1 m = \m) e 2 -*>-+^ \^ s ) 

and therefore the coefficient of \m,tfj s ) changes to 



(T) 



-e 



jeJ T 



E 



(T) 27ri(jii y 3i + ... + (j" 3 +p T+ i) ¥ > 3 + ...j 2 tip 2 t) 



cr -e 



But this is a trigonometric polynomial of powers Jr+i 
(recall eqn. (6)), and we can write it as 



^ 1} (fl= E CP* 



4 



if we define the coefficients a^ T+1 ) properly: 

m,s,j * * J 

5 (T+1) = a (T) 

m,s,(ji,...,j s ,...,j 2 t) m,s,(ji,...,j e -p T +i,---,j 2 t) 

for pt+i < js and otherwise. 

Finally we define Sm^i'fi) ■= Sm} s (<p) for the states 
for which the control bit is not set (m; t+1 = 0) and we 
can write 



W££U T ■ ■ ■ W£U \1,V>) = EE SZt^m K V< s > ■ 

m—0 s—1 

Now we use that the transformation Ut+i and the 
eigenvectors \ip s ) are fixed for all algorithms we consider 

(T I 1) 

and define w^ s;n t = (m,ip s \ Ut+i \n,ip t )- Then 



(/ T+1 ^C +1) ^)K^> 

= E S£ +1) (fl«£S, t |m,V>,>, (11) 



n.t 



which gives the following coefficient for \m^ip s ): 



E E 2 



(T+l) p 2^l<p (T+l) 



- e 



E E«;^ 5(T+1) 



(12) 



j e j7t+ i 

= : E a^e 2 -^:^^^). 



This completes the induction and establishes equations 
(9) and (10). 

Using the same argumentation as in equations (11) and 
(12) we can finally rewrite the state in equation (9) in the 
standard basis |fc) G {|0) , |1) , . . . , |2 C +* - 1)} through 

a™=Y(k\m,ip s )a( T+1 ) 

k,j t-^i * ' m,s,j 



and sj^P(ip) is of the same powers as Sm,s(<P), which is 
of powers Jj>. This proves equation (7) and (8). □ 

We now focus on the specific problem of phase esti- 
mation. The next lemma provides us with a necessary 
condition on the powers p\ , P2 , ■ ■ ■ such that a quantum 
algorithm with power queries can solve the phase estima- 
tion problem with precision e. 

Lemma 3. Any quantum algorithm estimating the phase 
ip of an eigenvector \q) of matrices Q from the class 

Q\q),t — {Q '■ Q * s a unitary t qubit transform, 

\q) is an eigenvector ofQj. 



up to precision e has to use power queries , , . . ., 
Wf T such that the set 

JT 

M T = {i-i'\i,i'= $> fc |tfc{i,...,T}} (13) 

keK 

has more than \M.t\ > elements. 

Note that this is both a condition on p\ , . . . , pr as well 
as on T: by cleverly choosing the powers pi, . . . ,pr we 
can get away with a smaller T, e.g. for pj = 2 3_1 we get 

M T = {-2 T + 1, -2 T + 2, . . . , 2 T - 1} , 

while for the choice pj = 1 we only get 

M T = {-T,-T+1,...,T-1,T}. 

Proof of lemma 3: For simplicity assume that 



2e 



— for some N e N. 
N 



(14) 



We will analyze the behavior of all possible algorithms 
for the phase estimation problem on a special subset of 
Q\q) t t- Fi x some arbitrary vectors \ip 2 ) , ■ ■ ■ , \ip2*) such 
that \q) , 1^2) , • • • , \ip2*) form an orthonormal basis and 
consider the following input: 



Q r :=e 2 ^\q) <g| + £e 



2iri<p s 



\i> s ) (ip s 



(15) 



s=2 



The phase ip we are interested in is ip = 2re for input Q r . 

Since the difference between the phases of the matrices 
Q r is 2e and we require e correctness, a measurement will 
yield states \k) from the distinct sets B r : 



B r = {k: \\2re - k\\ < 



!<«} 



(16) 



Depending on the number of qubits we use in our quan- 
tum algorithm, the sets B r can contain one or more 
states. 

By Theorem 2 we know that we can write the coeffi- 
cient of a state |fc) G B r after T queries as 



E 



a 



Ui,-,j 2 t)ej T 



(T) 

^(ji>-,i 2 i)' 



2ni(j 1 tf>+j 2 ip2 + ---+j 2 t'P2t) 



In this proof we are only interested in the behavior for the 
Q r . Therefore we can drop the dependence omp 2 , . . . , ip 2 t 
and let 

iec T 

where Co = {0} and 

£t ={ii ■ Uu ■ ■ ■ , j2* ) e Jt) 

={ji,h+PT ■ (ji,---,j2*) e Jr-l} 

= {J2Pk\KC{l,...,T}} 

keK 
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and the coefficients (note that I is fixed on the right side) state from B r is correctly returned from ps r (<£>), 



3 (T) 

Pk,l 



E 

(i,j2,...,j 2 t)eJ T 



a fc,(/ j2 ,...,i 2t ) e 



JV-l 



The probability pB r {f) of measuring a state from the 
set B r defined in (16) of all e approximations to tp = 2re 
is now given by: 

p B M :=£ \S ( P(V)\ 2 

k£B r 



(T) R (T) e 2ni{l-l')y 



= E E E OS 

fcGB r l€C T I'GCt 

= E E ^S^'-^ 

with the set M. T given by 

M T = {I - I' 6 jCt} 
and the coefficient 

„(T) _ V- -Y (T) V V ^ T , } /3 (T) 



(17) 



Li'ec T fe es r 

l — l'=m 



For illustration recall figure 2, which shows exactly one 
of these probability functions ps r (<£>) and also their highly 
oscillatory behavior. In the case of the phase estimation 
algorithm, B r = {\r}} and the figure shows pb 2 { 1 p)- 

We apply the Discrete Inverse Fourier Transform to 
PB r {f) (evaluated at the points ip = n/N for n — 
0, . . . , N — 1) and get for the fc-th coefficient 



E^(^)e- " 



n=0 



, r . 27rifcr 

>\PB r (jj)e « 



JV-l 



E^(]v )e 



n=0 
n^r 



(20) 



JV-l 



n=0 



since the probability pB r (<p) has to obey pB r (r/N) > | 
(recall the definitions of ps r and Q r ). If we knew that 
for the second term in (20) 



JV-l 



(21) 



n=0 

we could establish that 

JV-l 
n=0 



N 



E 



„(T) 
'Ir.ra 



m=k mod AT 



> o. 



We will show that this property, while not necessarily 
always true, will be true for at least most of the PB r (<fi)- 
There are N different possible outcome sets Bq , . . . , 
Bn-i- We know that for any <p — n/N all mutually 
exclusive probabilities of measuring a state from B r for 
r = 0, . . . , N — 1 have to add up to at most 1: 

JV-l 

Y p Br ( n / N ) < 1 f or n = 0, . .., N- 1, 

Let R< C {0, . . . , N - 1} be the set of all r for which 
(21) holds and R- the set for which it does not. \R < \ 
has to be greater than 1 since we can split 



JV-l 



(T) 2m(m-k)n/N 



n=0 
JV-l 



; E E ^:™ e " 

ra=k mod AT 



since for m ^ k mod TV 



AT-1 

£• 

n=0 



,2ni(m-k)n/N 



^ _ e 2-rri(m-k)N/N 
I _ e 2iri(m-k)l/N 



0. 



(18) 



(19) 



We can bound (18) by separating the part where a 



JV-l 



JV-l JV-l 



^=£i>££^( 



n— n— r— 

into the following parts: 



JV-l JV-l JV-l 

$>flr(;£)+ E Ep^(J)+ E E^(^) 

> -TV + Ii2< I + - |i£ I = -N + - \R* I 
4 4 4 4 1 1 

and therefore \R^\ < ±jV or \R<\ > |JV. 

Thus there is an r G i? < for which eqn. (21) holds and 



< 



N 



(T) 
'Ir.ra 



m=k mod N 



< N 



m=k mod AT 



(T) 
'/r,m 



(22) 
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for allfc = 0,...,JV-l. 

(T) 

This means at least N of the r]r,m have to be nonzero 
and thus PB r (<fi) from cqn. (17) must have at least N 
nonzero terms. In other words 

\M T \>N=^ (23) 

where we used the definition of N in (14). □ 

Lemma 3 now allows us to give a lower bound for the 
phase estimation problem. The numbers of queries T 
that is needed by any quantum algorithm is 0(log \). 

Proof of Theorem 1: From lemma 3 we know that the set 
Mt has to have \Mt\ > ^ elements. We can easily 
derive an upper bound on |A4t| ■ There are at most 2 T 
elements in the set 

= { £> fc |jrc{i,...,T}} 

keK 

and therefore we have \M T \ < \C T \ 2 < (2 T ) 2 = 2 2T . 
Combining our estimates for \Mt\ we get 

2 2T > \M T \ > 1 
and the number of queries T must grow like 

□ 



V. CONCLUSIONS AND EXTENSIONS 

In this paper we have obtained lower bounds for quan- 
tum algorithms that approximate the phase estimation 
problem through the new lower bound proof technique 
of frequency analysis. These lower bounds match the 
known upper bounds for phase estimation. 

The frequency analysis method can be used to give 
lower bounds for the Sturm-Liouvillc eigenvalue problem 
[7]. The results also extend to other forms of quantum 
queries, e.g. the query 

Q f \x) \y) = \x) \y © f(x)) = \x) \y + f(x) mod 2 m ) . 

An application of the frequency analysis method to these 
problems will be the subject of future work. 
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